TRI-PP-07-25 



Bridging over p-wave 7r-production and weak processes in 
few-nucleon systems with chiral perturbation theory 

Satoshi X. Nakamura 1 '!*) 

1 Theory Group, TRIUMF, 4 004 Wesbrook Mall, Vancouver, BC V6T 2A3, Canada 

(Dated: May 14, 2008) 

Abstract 

I study an aspect of chiral perturbation theory (%PT) which enables one to "bridge" different 
reactions. That is, an operator fixed in one of the reactions can then be used to predict the 
other. For this purpose, I calculate the partial wave amplitude for the p-wave pion production 
(pp — > pnir + ) using the pion production operator from the lowest and the next nonvanishing orders. 
The operator includes a contact operator whose coupling has been fixed using a matrix element of 
a low-energy weak process (pp — > de + u e ). I find that this operator does not reproduce the partial 
wave amplitude extracted from experimental data, showing that the bridging over the reactions 
with significantly different kinematics is not necessarily successful. I study the dependence of the 
amplitude on the various inputs such as the NN potential, the 7T-/VA coupling, and the cutoff. I 
argue the importance of a higher order calculation. In order to gain an insight into a higher order 
calculation, I add a higher order counter term to the operator used above, and fit the couplings to 
both the low-energy weak process and the pion production. The energy dependence of the partial 
wave amplitude for the pion production is described by the operator consistently with the data. 
However, I find a result which tells us to be careful about the convergence of the chiral expansion 
for the pp — > pmr + reaction. 
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I. INTRODUCTION 



Applications of the chiral perturbation theory (%PT) to few-nucleon systems have been 
of great interest in the last decade. For reviews, see, eu)., Refs. [lj. Particularly, a compre- 
hensive review for the tt production is given in Ref. [2J, and for the electroweak processes 
in Ref. [jj. A powerful aspect of xPT is that it provides us with a bridge between different 
reactions in a model-independent way. This means that coupling constants (the so-called 
low-energy constants, LECs) fixed using experimental data for one of the reactions can then 
be used in the other. An interesting interaction in this context is 

C = -2dxN ] S ■ uNN^N , (1) 

with 

/> M = -r a d^TT a - e 3ab V^n a T b + f n A^ H . (2) 

The spin operator is S, and the external vector (axial) current is (A^). The constant d\ 
is a LEC. This contact interaction between the nucleons contributes to the three-nucleon 
forcejj], the p-wave 7r-production (pp — ► pn7r + [5j, pp — > dn + ), the radiative pion capture 
on the deuteron (n~d — > 'ynn)^, Q], and the weak processes in few-nucleon systems, such 
as tritium /5-decay and the proton-proton fusion (pp — > de + i> e )\^, the neutrino-induced 
disintegration of the deuteron [y e d — > e~pp, vd —> upn)^, and the muon capture by the 



deuteron (fi~d — > v^nri) [10|] . The contact term in Eq. ([T]) makes the connection among 
these reactions, which may be referred to as the two-body analog of the Goldberger-Treiman 
relation, as stated in Ref. [6|. If all couplings except for d 1 have been fixed using experimental 
data for ttN and iViV elastic scattering, then one can fix d\ using one of the above reactions, 
and predict the others. 

There have been several such calculations which I will refer to as the "bridging program" . 
One of them was done by Park et al. where they fixed d\ using the experimental tritium 
/3-decay rate, and calculated, with no free parameters, the weak proton capture by a proton 
(or 3 He). In this case, all the reactions are low-energy weak processes, and the kinematics 
are relatively similar. Therefore, the bridging program is expected to work well. In another 
work due to Hanhart et al. [5|], the authors calculated the partial wave amplitude for the 
p-wave 7r-production {pp — > pnn + ), and showed that the use of d\ fixed by three-nucleon 
observables ^| also consistently reproduces the partial wave amplitude extracted from ex- 
perimental data Although this result seems to be satisfactory, the bridging program 
in this work was nol dee fu,.y consistently. This is because they used a uuelea. force 
which is different from the nuclear force used for fixing d\ in Ref. 0] . Because of the short- 
range nature of the d\ term, the d\ value is largely dependent on the choice of the nuclear 
force. Hanhart et al. also showed that the p-wave 7r-production amplitude is rather largely 
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dependent on the contribution from the d± term. This indicates the importance of careful 
treatment of the d\ term in the calculation; the nuclear force and the d\ value have to be 
consistent. Finally, I mention the bridging program done by Gardestig et al. [6], where they 
fixed the d\ value using the matrix element of pp — *> de + u e and then used it in calculating 
observables for n~d — > 77m. The two reactions are similar in kinematics of the NN sector, 
but are rather different in the momentum transfer, which might have to be taken care of. 
They are interested in extracting the neutron-neutron scattering length from the n~d — > 'ynn 
reaction and, for that purpose, they are interested in the shape of the spectrum rather than 
the absolute value. They found that the use of the d\ term fixed in the above manner signif- 
icantly reduces the dependence of the shape on R which is the matching point between the 
long range one-pion-exchange potential and the short range square well potential. Although 
this result is a success of the bridging program, it is still interesting to study the absolute 
value of the cross section in order to more rigorously test the power of xPT. 

In this work, I would like to more seriously investigate how reliably the bridging program, 
an important aspect of %PT, works. I believe that my investigation is important because 
there has been sometimes an argument which supposes that the bridging program works, 
even though it has not been quantitatively confirmed that a bridging program over reactions 
with considerably different kinematics works. For this purpose, I calculate the partial wave 
amplitude for the p-wave 7r-production in NN collision {pp —>■ pmr + ), with d\ fixed by 
an observable of a low-energy weak process. This obviously provides a stringent test of 
%PT, because the two reactions are strong and weak processes, and are quite different in 
kinematics. My program is as follows. I extend the operator of Ref. [§J by including the A, 
which is known to play an important role in the p-wave 7r-production, and re-fit d\ to an 
observable of a weak process. Here, I fit d\ to a matrix element of the axial vector current 
used in Ref. jsj: the kinematics is that for the pp — > de + v e reaction. The 7r-production 
operators I use are the same as those used in Ref. |5| , except that the A is treated differently. 
In Ref. |^| . a coupled-channel equation was solved and the A is included in the nuclear wave 



functions 12|. I use nuclear wave functions including only the nucleons as dynamical degrees 
of freedom. I take account of the A by including it in the 7r-production operator. 1 The 
7r-production operator I use is based on a counting rule proposed in Ref. p, [3], which is 
different from Weinberg's counting[14], and the large initial on-shell momentum [~ ^/m^m vi 
m N { m ir) '■ the nucleon (pion) mass] is considered as a characteristic scale of the system. I 
consider the operators up to next-to-leading order [NLO, ©(m^/rajy)]. 2 



1 Thus I treat the A perturbatively and do not fully take account of some non-perturbative effect of the 
A in wave functions. A calculation with a fuller account of the A might be worthwhile doing to see a 
difference. 

2 The operator used in Ref. [8[ is based on Weinberg's counting while the operator used in Ref. [5[ is based 
on the counting proposed in Ref. [3, Therefore, one may wonder whether the two operators contain 
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It is very interesting to see whether the operator constructed using the low-energy weak 
process can reasonably describe the 7r production. I will use several combinations of the 
inputs (the tcNA coupling, the NN potential, cutoff) in my calculation. Even though the 
di value is fixed for each combination of the inputs so that the low-energy weak process is 
reproduced, one may expect a dependence of the tt production amplitude on the combination 
because of the rather different kinematics between the two reactions. I will study such a 
dependence. In fact, we will see that this bridging program is not successful when working 
with the NLO operators; the partial wave amplitude of the 7r production based on x?T is not 
consistent with the data. I will argue the importance of going to a higher order calculation. 
In order to, even roughly, explore a result of a higher order calculation, I add a higher order 
counter term and study a consequence. 

This paper is organized as follows. In Sec. [Til 1 present the chiral operator for the p-wave 
7r-production up to the order I work with. I discuss the determination of the d± value using 
the low-energy weak process in Sec. IIHI In Sec. IIV} I perform a multipole expansion of the 
operator, and express the cross section in terms of the partial wave amplitudes. Then, in 
Sec. El I present my result for the p-wave tt production amplitudes and compare them with 
the data. The result obtained with the higher order counter term is also presented. Finally 
I summarize this work in Sec. IVII 



II. CHIRAL OPERATOR FOR p-WAVE vr-PRODUCTION 

I present expressions of chiral operators which contribute to the p-wave 7r-pro duct ion. 
I basically follow Ref. [5J, except for the treatment of the A degree of freedom, and the 
way how the high momentum components are cut off. I start with the chiral interaction 



Lagrangian given in Ref. 15J. By keeping terms relevant to my calculation, I have 



42 = -T^T • (W >< *)N + |f N] (™ ■ V7T)iV 

l>A N\TS- V7r)A + h.c. , (3) 



and 



V, 



42 = ir^-n^T ■ (tt x Vtt) ■ ViV - %N\Vtz) 2 N 
-N ] -^-a ■ Vtt x Vtt • tN — iVV • ttg? • VN 



different mechanisms and thus the d\ value is different in each case. I will argue in Sec. IIIII that the d\ 
value should be the same in the both cases, up to the order 1 am working. 
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(a) 



(b) 



(c) 



(d) 



(e) 



FIG. 1: The p-wave 7r-production operators up to NLO in xPT. Dashed lines denote pions and 
the solid lines nucleons, the double lines A. Vertices without (with) the shaded box arise from the 
leading (subleading) order Lagrangian. 



d 2 



2/. 



Vtt x N ] arN ■ N^arN + h.c, 



(4) 



with C4 = C4 + 4^—- The nucleon's spin (isospin) operator is a (t), while the NA transition 



spin (isospin) operator is S (T). I use the pion decay constant /„. = 93 MeV and the 
nucleon mass, = 939 MeV. Regarding the parameters (<7A) ^M.j c 3> c 4)> I follow Ref. [161 ]. 



where Krebs et al. constructed a chiral nuclear force including the A explicitly. The axial 
coupling constant is qa 1.27). For the 7riVA coupling, I use two choices: Ha.— 2.10 



from the A-decay width 171]. and Ha = 3gA/V% = 2.68 from large No- 3 In Ref. 16j, the 
authors calculated the s- and p-wave irN scattering threshold parameters at next-to-leading 
order with taking account of the A, and fit the couplings (03, C4 and others) to the analysis 
of Matsinos[18]. The result is c 3 = -1.87 GeV" 1 , c 4 = 1.87 GeV" 1 for h A = 2.10, and 
C3 = —0.79 GeV -1 , C4 = 1.33 GeV -1 for h-A = 2.68. The remaining unknown LECs are d\ 
and d>2 which will be determined in the next section. 

I use the following 7r-production operator derived from the above interaction Lagrangian. 
The leading order [LO, 0(1)] operator is the one-body direct production of the pion with 
the isospin state a off the nucleon [Fig. [0(a)], 



O 



IB 



. 9 A 



(2n) 3 S^(p' 2 - P2 )t?ct 1 q n + (1 <- 2) 



(5) 



where q n is the momentum of the emitted pion. The quantity Pi (pQ is the momentum of 
the incoming (outgoing) i-th nucleon. Another LO mechanism I consider is the A-excitation 
followed by the n emission. In Ref. j^J, the authors used the wave function which explicitly 
includes the A component, and considered the one-body operator which produces the pion 
with the A deexcited to the nucleon. Because I use nuclear wave functions with only 



3 The definition of the 7T7VA coupling constant here is different from that of Ref. [16] by a factor of 2: 
h A (this work) = 2h A (Ref. 
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the nucleonic degrees of freedom, I alternatively use a two-body operator in which the A is 
excited either by the 7r-exchange [Fig.QJb)], or by a contact interaction between the nucleons 
[Fig. (He)]. The former is given by 

Q j_ 9Ah\ ct> ■ k (4r 2 a fc - (n x T 2 ) a cr l x k ) ■ q n 

An ' 36 /3 m» + k 2 mA -m N -pl/m N + ( P ' + qj2y/2 f i + 1 j ' U 

where (= 1232 MeV) is the mass of the A, and p Q is the initial on-shell relative NN 
momentum. The quantities, /i, k, p' are respectively defined by \i = vn^vn^lira^ + wa), 
k = p 2 — p 2 and p' = — p' 2 )/2. I assume the equal energy sharing between the nucleons 
and = (3m 2 — g^)/4, m„- = 138 MeV. In Appendix, I explain how I treat the energy 
denominator in Eq. ([6]) in my calculation. An expression for the diagram in Fig. []Jc) is 

i h A C 2 A (4 r 2 a (cr 2 x k) x fc., - (n x r 2 ) a <Xi x [(cr 2 xfe)x fe]) ■ q w 



O 



ACT 9 / w m A - m N - p 2 /m N + (p' + q 7T /2) 2 /2fi 

+ (1 - 2) , (7) 

where I have used the contact NN — ► iVA interaction with two derivatives, 

£>nn—>na = — «C 2 VA A t T ■ (V + V)5 - 5a • (V + V)) iV- iVt (?x(V + V)) riV .(8) 

It is noted that a contact NN — > iVA interaction without derivative does not contribute to 
the transition under consideration. I use the C 2 A value taken from the resonance saturation 



of the p-exchange NN — ► NA potential used in several phenomenological models 19] : 



h ^ NA = 0*1 ^ + *>? k • < 9 > 

where the pNN vector coupling is g P NN {.Q 2 p n n I ^ = 0-5) and the tensor coupling is k p 
(= 6.6). The factor at the end, 1/m 2 (m p = 770 MeV), is the static limit of the p-meson 
propagator. Although Eq. (J5J) is not a general contact NN — > NA interaction with two 
derivatives, one may take account of the most important part of the NN — > NA contact in- 
teraction by invoking a meson-exchange model, and the remaining part may be phenomeno- 
logically absorbed by the contact d{ (i = 1, 2) terms. Obviously, my treatment of the A 
is rather phenomenological. However, a construction of a nuclear force from a chiral La- 
grangian with the tt, N and A, which is yet to be done 4 , is necessary to determine the 
C 2 A value. I believe that my treatment is one of what one can do best for the moment, and 
expect a fully consistent calculation in future. 

Next I discuss next-to-leading order [NLO, 0(m n /mjv)] terms which consist of four types. 
One of them is the recoil correction to the LO terms. The recoil corrections to the one-body 
term [Eq. (J5J)], 7rA term [Eq. (jB])], and contact- A term [Eq. are respectively given by 

1B , reC oii = -t-^ r (27r) 3 5^(p' 2 - p 2 )r 1 V 1 • ( Pl + p[) + (1 ~ 2) , (10) 



Such a nuclear force for the peripheral wave has been constructed recently [1 
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Oatt. 



recoil 



X 



w w g A h\ cr 2 • k 



72 m w /3 



m 



/2 



+ fc 2 



4r 2 a fc • (pi + pi) - (n x T 2 ) a [(<r a x k) ■ (p a + pi) + 2iA; 2 



m A -m N - pl/m N + (p' - q^) 2 j2\i 



+ (1 



ACT.recoil 



1 



18 m A -m N - p 2 /m N + (p' - q 7r ) 2 /2/j, 

{r 2 a [-2A; 2 fc ■ (<n x cr 2 ) + 4((<r 2 x fe) x k) ■ (pi + p[) 
(n x r 2 ) a (o-i x ((<x 2 x fc) x fe)) ■ (px + pi)} + (1 <-> 2) 



;i2) 



where ^ is the energy of the emitted pion. The second type of the NLO operator is a pion 
rescattering through either the C3 term, or the C4 term or the Galilean correction to the 
Weinberg- Tomozawa term [the first term in Eq. (jlj)]. They are graphically represented by 
Fig. [D(d), and their expressions are given by 



a 



■ Cz9A 



(To ■ k 



fl m' 2 + k 2 



k ■ + (1 <-> 2) 



(13) 



°« = -^^Tk2^ >< fc ) ■ ftrfa x r 2 f + (1 - 2; 



(14) 



WT, Galilean 



(T 2 • 



C« w " *0 • (Pi x r 2 ) a + (1 ~ 2) . (15) 



16m n fl m' 2 + k 2 

For a convenience, I decompose the last term as OwT,Gaiiiean = Owti(g) + Owt2(g) with 5 

9 a cr 2 - k 



O 



WT1(G) 



o 



WT2(G) 



16m N f% m' 2 + k 2 

9 a <r 2 k 
16m n fn m n + ^ 2 



(2p' + k)-q 7V (r 1 XT 2 ) a + (l<->2), (16) 
(-k 2 - 2p' ■ fe)(n x r 2 ) a + (1 <-> 2) . (17) 



5 When convoluted with the wave functions, Owt, Galilean gives a non-vanishing amplitude in the soft pion 
limit. (The tree amplitude vanishes in the soft pion limit.) This is not consistent with chiral symmetry. 
This problem is similar to that found in Ref. [201 ] in the context of the s-wave pion production in the 
NN collision. The solution to this problem was given in Ref. 21]. Probably, the problem here is also 
resolved by the same mechanism found in Ref. [21(; e.g., one loop diagram formed by Owt, Galilean and the 
one-pion-exchange potential is (partly) canceled out by other irreducible pion loop diagrams, leaving a 
contribution consistent with the chiral symmetry. One may take some prescription to maintain the chiral 
symmetry. However, I use Eq. (|15p without modification because this term gives a small contribution (a 
few percents) to an amplitude for the pp — > pmr + reaction; the modification will not significantly change 
results. 
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The third type of the NLO terms is a pion rescattering via the Weinberg- Tomozawa term: 

°« = - x ^ + (1 ~ 2) ' (18) 

Finally, the fourth type is a pion emission from the contact terms [the d\ and d 2 terms in 
Eq. ®, Fig. m^e)]: 

O d = -i {j^i + x °"2)(r 1 x r 2 )^ .q„+ (1 <- 2) . (19) 

Starting with the operators presented above, I take the following procedure to obtain 
the momentum space (p-space) operator contributing to the ^So ^ 3 S\- 3 D 1 transition. At 
first, I perform Fourier transformation of the above operators to represent them in the 
coordinate space (r-space). In the Fourier transformation, I multiply a Gaussian cutoff 
function, exp (— A; 2 /A|>), to the operators other than the one-body operators; I use A G = 
2 GeV. In the r-space, I perform the multipole expansion of the operators, and evaluate 
the matrix elements for the spin, isospin and angular parts of the operators. The remaining 
radial part of the operators is converted to the p-space using Fourier transformation, and 
I obtain the radial part of the p-space operator, 0(p',p), where p (p') is off-shell relative 
momentum of the incoming (outgoing) two nucleons. I classify the operators into two groups: 

Group I : Oatt, Oact, C3 , C4 , O d 

Group II : Oib, OiB,recoih ^A^recoU) OACT.recoil , OwT.Galilean, 0\yT 

I introduce a sharp cutoff A for the operators belonging to Group I such that 

0\(p',p) = 0(p',p) , for p < A and p < A , 

0\{p',p) = , otherwise . (20) 

I do not apply the sharp cutoff to the operators of Group II. As we will see in the next section 
where the LECs d{ are determined using low-energy weak processes, the renormalization of 
di only takes care of high momentum components (larger than A) of the operators in Group 
I. The LECs di contain the same physics for both the reactions bridged. Therefore, I retain 
the high momentum components of the operators in Group II in my calculation. 6 As 
representatives, I choose A = 500, 600 and 800 MeV. I use the p-space operator explained 
above because of its usefulness for the renormalization group (RG) analysis which will be 
done later. 



6 One may also apply the sharp cutoff to the Group II operators, which amounts to cutting off a higher 
order effect. 
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For the purpose of the multipole expansion of the ir production operator, which will be 
done in Sec. HV] I express the above n production operators as follows: 



O x = -Ox-Q, X= IB, Att, ACT, c 3 , c 4 , WT1(G), d (21) 
Y = Y Q°, Y = recoil terms, WT2(G), WT (22) 

where Q = iq^llf^ and Q° = iuj 7r /2f n . I will refer to Ox as the spatial component of the ir 
production operator while Y as the time component. 

III. DETERMINATION OF d t 

In this section, I determine the LECs di in Eq. ffT9l using an observable of a low-energy 
weak process. The di terms contribute to the 1 5 — > 3 5*i transition and appear as a single 
linear combination, d = d± + 2d2- Because di and g?2 are not separately determined by 
considering the 1 So ^ 3 Si transition only, I determine d rather than d\ and c?2 individually. 

In Ref. [8|], the authors fixed d using the experimental data of the tritium /3-decay rate. 7 
They derived the axial vector current contributing the tritium /3-decay from the chiral La- 
grangian. They did not explicitly consider the A as a degree of freedom. The spatial 
component of the axial vector currents used in Ref. [s| are obtained from the one-body op- 
erator [Eq. (j3j)] and Group I operators without the A [Eqs. (|T3|) . (THj) and (fl9l) ]. with the 
factor of {iqn/^fn) being eliminated and m' n replaced by m v . Accordingly, the parameters 
(cj) used in Ref. [8j are different from those presented in the previous section. Although 
they additionally included some other operators which give a negligible contribution, I do not 
consider them in the following. They multiply a Gaussian cutoff function, exp 
to the operators, with A x = 500, 600 and 800 MeV. For each choice of A x , they fixed d 



22 



so that the tritium /3-decay rate is reproduced. They used the AV18 iVTV potential 
supplemented by the Urbana-IX three-nucleon potential [23) when calculating the matrix 
element for the tritium /9-decay. 

Here, I need to re-fit the value of d in "my" axial current operator. My operator is the axial 
current operator used in Ref. [s| plus the A-excitation current. The A-excitation current is 

obtained from Eqs. ([6]) and ([7]) by eliminating the factor of (iq n /2f n ) and m' n — > m n . The 

f L 

way of cutting off the high momentum states is also different from that used in Ref. [8j, as 
has been discussed in the previous section. I do not directly use the tritium /3-decay rate 
but take an indirect way to fix d, as discussed in the next paragraph. 

I start with a benchmark calculation to which the d value in my operator is fitted. For that 
purpose, I use the same spatial axial current operator used in Ref. ^ (the same couplings 



The authors fixed dimensionless constant cLr rather than d. The two quantities are related through 
d R = (m N f%/g A )d + c 3 /3 + 2c 4 /3 + 1/6, with c 3A = c SA /m N . 
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TABLE I: Dimensionless contact coupling, d(= m^f^ d). The first column is the sharp cutoff 

value. The first row specifies the NN potentials used. For each NN potential, the left side is d 
for Ha = 2.10, while the right side is for Ha = 2.68. 

A CD-Bonn AV18 Reid93 Nij I N 3 LO 

500 0.18 0.06 -0.01 -0.17 -0.48 -0.64 0.29 0.15 -0.03 -0.20 

600 0.62 0.42 0.77 0.47 0.17 -0.13 0.93 0.67 0.60 0.32 

800 1.74 1.32 4.36 3.32 3.16 2.11 2.75 2.13 1.18 0.81 



and the same cutoff but without the negligible operators) to calculate a matrix element for 
the 1 S'o ^ 3 Si- 3 Di transition in two-nucleon system. This matrix element is the benchmark 
to which d is re-fitted so that "my" operator reproduces the same matrix element. In the 
benchmark calculation, I use the AV18 potential for consistency, and choose a kinematics 
with T l NN = 0.5 MeV (T^ N : the initial on-shell kinetic energy of the relative motion), the 
deuteron final state, and q = 2.5 MeV (q : the momentum transfer from the two-nucleon 
system to the external current). This kinematics is for the low-energy pp — > de + u e reaction. 
I use the proton-neutron interaction to generate the initial wave functions, and therefore 
do not consider the Coulomb interaction and the isospin violation effect. I calculate the 
matrix element with different choices of the Gaussian cutoff, A x = 500, 600 and 800 MeV. 
I average the three matrix elements which have a slight cutoff dependence, and regard the 
average as the benchmark. When I calculate the matrix element of my operator, I use several 
combinations of the NN potential, the 7riVA coupling (Iia) and the corresponding cs values, 
and the^sharp cutoff value (A). I use the following NN potentials: the CD-Bonn 24J, the 



AV18[22|, the Reid93[2a|, the Nijmegen I[25j and the chiral N 3 LO[26| potentials. 8 For each 
combination of these inputs, the d value is fitted to the benchmark, and the result is given 
in Table [U In the table, I show the dimensionless coupling d{= m^f^ d). Although the d 
value is adjusted so that the benchmark is reproduced irrespective of changing the inputs, 
it would be expected that the 7r-production amplitude, evaluated in a significantly different 
kinematics, have the dependence on them. I will study the dependence later. 

Before closing this section, I discuss the issue mentioned in footnote [2J In order for 
the bridging program to be meaningful, the d term has to (implicitly) include the same 
mechanisms (up to the external probe) for both of the reactions bridged. If the d term for one 
of the reactions includes a mechanism which is explicitly considered in the other, the d value 



I take the hybrid approach where a matrix element of the chiral operator is sandwiched by wave functions 
obtained with a phcnomcnological nuclear force. This approach is also referred to as more effective effective 



field theory (MEEFT). An argument for employing MEEFT is given in Ref. [27 1. 
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should be different in each case. In the case under consideration, the axial current operator 
and the pion production operator are based on different counting schemes, and therefore 
there is a concern that the two operators explicitly contain different mechanisms. In the 
following, I will argue that the d term includes the same physics for both of the operators, 
up to the order I am working. As we have seen in the above, all of the mechanisms for the 
axial vector current are included in the pion production operator (up to the external probe). 
However, the pion production operator contains some mechanisms which are not included 
in the axial vector current; the recoil corrections, WT and WT(G). Among them, WT and 
WT(G) do not contribute to the weak process, and thus the d term fixed in the weak process 
does not contain these mechanisms. The recoil correction of the one-body axial current [the 
counterpart of Eq. (1101) ] has been considered in Ref. [sj], which means that the d term fixed 
in Ref. does not contain this mechanism. Although I did not consider this mechanism 
when fixing d, I am safe because this mechanism gives only a negligible contribution to the 
matrix element considered. The recoil corrections of the A axial current [the counterpart 
of Eqs. ffTTj) and ffT2l) ] have not been considered in Ref. However, the d value does not 
depend on the inclusion of these mechanisms because they give a negligible contribution to 
the matrix element. (In fact, these recoil corrections should be captured by another higher 
order counter term.) Therefore, up to the order I am working, the d term contains the same 
physics for both of the reactions bridged, and I can use the d term fixed in this section for 
calculating the pion production amplitude. 

IV. PARTIAL WAVE AMPLITUDE AND CROSS SECTION 

In this section, I will express the cross section for the pp —>■ pmr + reaction in terms of 
partial wave amplitudes. 9 For this purpose, I perform the standard multipole expansion 



of the it production operator 28j], and a partial wave expansion of the initial and final NN 
scattering wave functions. Then, I use the partial wave amplitudes to express the transition 
amplitude, and subsequently the cross section. I perform these expansions in the r-space, 
which will be followed by the conversion into the p-space. 



A. Partial wave amplitude 

The multipole operator for the time component of a it production operator is defined by 

T£ M (0) = J dxjj(qx)Y JM (x)O {x), (23) 



I will work with the center-of-mass system of the initial pp system throughout this work. 
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where q = —q w (q = \q\), jj(qx) is the spherical Bessel function of order J, and x = 
x/\x\. An r-space operator is O (x), the dependence on the center-of-mass coordinate 
being eliminated; e.g., 0°(x) = O°(x)\f{^ . For the spatial component of the tc production 
operator, I show only the longitudinal multipole operator because the electric and magnetic 
multipole operators give vanishing contribution for the case in question: 



T J L M {0) = - dxV[jj(qx)Y JM (x)]-0(x). 
q J 

The transition amplitude for the pp —>■ pmc + reaction is written by 



J2<d° y (x)q -j:<d x (x).q 

Y X 



(24) 



\1>i) 



(25) 



where ipi (if>f) is the nuclear wave function for the initial (final) state, whose r-space repre- 
sentation will be given later. In the summations in Eq. (j25j) . X (Y) takes various components 
of the 7r production operator specified in Eq. (1211) [Eq. fl22l ]. The four- vector Q has been 
defined in Eqs. fT2T|) and (1221) . Using the multipole operators presented above, I can rewrite 
Tji as 



T fi = E 47^-1)^/1 

Jo Mo 



T, 



Jo M r\J —M 



C 



'Q 



c 



JoM r)Jo—M 



+ T J L ° M °Q 



(26) 



with 



n JM 



Qt M 



Y JM (q) Q°, 



J 



2J+1 



■ J-1JM 



(9) 



/ J+l 
2J + 1 



Yj +lJM {q)\ Q 



(27) 
(28) 



and Yjlm{q) are the vector spherical harmonics. 

Now I proceed to the partial wave expansion of the iViV wave function. An iViV scattering 
wave function with the relative momentum p, the third component of the spin (isospin) of 
i-th nucleon being Sj (rj) is expanded as follows: 

V>(r) = E 4vr(l/2, Sl , 1/2, s 2 \Sfi)(l/2, r 1; 1/2, r 2 |T, T 3 )(LmSfi\ JM)^ L F L * m (p)^(r), (29) 



where the index a collectively denotes the quantum numbers of a partial wave; a = 
{J, L, S,T} where J,L,S are the total, orbital, total spin angular momenta of the iV\/V 
system, respectively and T is the total isospin. The partial wave function (ip a ) is expressed 
as 

_l\L+S+T 



1 



r- J2yL'Sj(r) Ri LS (r) T] TiT 3 
V2 v 



y LS j{r) = \Y L (r)®xs] 



(J) 



(30) 
(31) 
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where the two-nucleon spin (isospin) wave function is denoted by xs (Vt)- The radial part 
of the above wave function is normalized, in the plane wave limit, to be 

Ri, L; s(r)^jL(pr)5 L>L , (32) 

With the multipole operators and the partial waves presented above, I express the tran- 
sition amplitude in terms of the partial wave amplitudes: 

Tfi = e E E (-i) Mo ^ Jo+Li - Lf ^=yum Y Lf , mf (p f ) 

ai,mi<xf,mf J ,M V / ^ 



x (1/2, si,*, 1/2, 8 2 ,i\SifH) (1/2, r M , 1/2, T 2ji \T h if) (L.m^l J^M,) 

x (1/2, s ltf , 1/2, 52,/!^) (1/2, ri f/ , 1/2, r^T,, T|) {Lfm s S f iif\JfM f ) 

x^M^MolJfMj) £ (T a J °> Q a J °" M ° , (33) 

a=C,L 

with the suffix z (/) indicates the initial (final) state. In the above equation, I used the 
abbreviation 

<r/-> = ^ Q/ ||T/°||^ S ) , (34) 
for the reduced matrix element defined by 

(J^MflT^lJ^Mi) = -=L (JMJoMolJfMf) (J f \\T^\\J t ) . (35) 

y 2J f + 1 

B. Cross section 

The unpolarized cross section for the pp — > prm + reaction is 

^ = E — TTls^-^^ " ?t - Pf)\T f i\ 2 dq v d Pl jdp 2 j , (36) 

where Ya / indicates the average (summation) over the initial (final) spin and isospin states 
of the two nucleons. The quantities P, and Pf are the initial and final four total momentum 
of the two-nucleon system; P 2 = (P ) 2 — P 2 . The relative velocity between the initial two 
nucleons is denoted by v re [. 

I derive the pion angular distribution, retaining only partial wave transition amplitudes of 
interest. I am primarily interested in the — > 3 Si transition amplitude, (ip3 Sl | |T/°| \i^is )i 
where the d term plays an important role. [It is noted that ip^Si contains the 3 Pi component, 
as seen in Eq. fl30|) . and thus the 1 5'o — > 3 D 1 transition is also included in (ilj3 Sl \\T^ \\ipi So ) ■] 
However, I also retain the 1 D 2 — > 3 Si transition amplitude for a later convenience. I retain 



13 



multipole operators with rank one (J Q = 1) which dominantly induce the transition. I 
integrate over the final two nucleon momenta to obtain: 



lcT = dq * — h T2 — l M l ' ( 3T 

all^ Jo 1d7t p io n j£ 



with 



\M\ 2 = \ £ K^JI^-g^ll^jr + Ps^os^QK^JlT^-g^ll^)! 2 

4 a i =i5o, 1 D2 V4 



(^ 35l | \T l c - q n T l L \\^ D2 ){^ Sl 1 17* - g w T2| |^o)*l J • (38) 



^Re 

The scattering angle of the pion with respect to the direction of the initial proton is denoted 
by 6 n , and Pzix) is the Legendre function of degree two. The initial one nucleon energy is 
E p and the maximum magnitude of the pion momentum is given by 

J {E p + m N + m 7T /2)(E p + m N - m n /2)(E p - m N + m n /2)(E p - m N - m n /2) 
Qn = 1 p • ( 39 ) 

Hyp 

I will compare my calculation with "partial wave amplitudes" extracted from experimental 
data in Ref. llj. However, the "partial wave amplitudes" are actually different from the 
ordinary one, (T^° — q n Tl°). Therefore, I establish the relation between them so that the 
comparison makes sense, which will be done in the next paragraph. 

In Ref. 11] . six lowest partial wave amplitudes are assumed to contribute to the 
pp — > pmr + reaction in the energy region under investigation, and are extracted from the 
data. In extracting the amplitudes, they treated the final iViV system as the "pseudo- 
deuteron" which is the iViV scattering state with the relative motion integrated over; the 
pseudo-deuteron may have the angular momentum different from the deuteron. This means 
that they parameterized the data using a formulae in which the pp — * pmr + reaction was 
regarded as the pp — > u d"ir + reaction; "<i" is the pseudo-deuteron. More specifically, they 
parameterized the pion angular distribution in the pp — > pmr + reaction using the formula: 

da 



dVt^ 



C + C 2 P 2 (cos6 n ) , (40) 



with 



_ |a | 2 |ax| 2 |a 2 | 2 |&o| 2 N 2 N 2 ( , 

|a 2 | 2 Re[a a;] 

where a (a 2 ) is the partial wave amplitude for 1 S — > 3 S\ ( 1 D 2 — ¥ 3 S\); see Table IV of 
Ref. for the other partial wave amplitudes. By comparing Eq. (j4"Ul) with Eq. fl37|) . I can 
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find the relation between the two differential cross section formula. To see the relation more 
clearly, I denote the theoretical counterpart to a (a 2 ) by a (a 2 ) and rewrite Eq. (I3TI) as 



da \a \ 2 \a 2 \ 2 ( \a 2 \ 2 Re[a a 2 }\ 

in; = — + — + {-4 ^-J Pa(008 ^ 



(43) 



with 



|«o| 2 



l fl 2| 2 



(^sjKc-^TillViSo)! , 



Re[aoa 2 ] 



T E p^l + P flp'£ u , ,. Tl Tl|l/ N|2 

aq-x — 

o 16vr p uj n f£ 



(44) 
(45) 



x Re 



(46) 



I distinctly used a and a because they are not necessarily the same. 10 Which (a or a) should 
be compared with ao and a 2 from the experimental data ? I take the following way to find 
a solution. 

At first, I factorize out the phase coming from the initial state interaction as 



a 



(T 



a 2 



S-2 7.1 



u 2 i 



(47) 



where 5 (5 2 ) is the phase shift of the { X D 2 ) partial wave scattering. I choose a' and 
2 to be real. The same factorization applies to ao and a 2 . Then, as Step (i), I calculate 
and \a' 2 \ from Eq. ([44]) and Eq. ([45]) . respectively. Next, as Step (ii), I solve a coupled 



a 1 



equation consisting of Eq. (j4"6"j) and |ao| 2 + |a 2 | 2 = |ao| 2 + |a 2 | 2 [the r.h.s. is from Step (i)], 
thereby finding a set of solutions, a and a 2 . Step (ii) is closer to the way how ao and a 2 are 
extracted from the data; each amplitude is extracted from the sum of the amplitudes (data). 
My numerical result shows that the two sets of the solutions obtained in Step (i) and (ii) are 
essentially the same in the absolute value. [Step (i) cannot fix the sign.] 11 Therefore, I will 
employ the solution from Step (ii), in which the relative phase between a and d 2 can also 
be fixed. (Hereafter, I do not distinguish between a and a, and denote them by a.) Still, 
the overall phase of a has not been fixed. In order to make the comparison with the data 



10 I use the symbol a (a) in a generic sense, referring to both ao and a 2 (5o and 02). 

11 Strictly speaking, a obtained in my analysis is not necessarily the same as a from Ref. In order 



to obtain a theoretically, one calculates the analyzing power and the pion angular distribution for the 
PP pmt + reaction, taking account of all non- negligible partial wave amplitudes. Then a is extracted 
from those observables in the same way as done in Ref. However, the fact that a obtained in Step 
(i) and (ii) are essentially the same would indicate that my result would not change drastically even if I 
took the "ideal" procedure to calculate a. 
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FIG. 2: The transition amplitude a 2 for pp — > 
pmr + . The chiral NLO it production operator 
and the CD-Bonn NN potential are used. The 
solid curve is obtained with Ha = 2.10 while the 
dashed one with Ha = 2.68. Experimental data 



are from Ref. 



FIG. 3: The transition amplitude a 2 for pp — > 
pmr + obtained with the chiral NLO ir produc- 
tion operator. I employ hj\ = 2.10. The 
solid, dashed, dotted and dash-dotted curves 
are respectively a 2 obtained with the CD-Bonn, 
AV18, Reid93, Nijmegen I NN potentials. The 
result with the N 3 LO potential almost exactly 
falls on the dash-dotted curve. Experimental 
data are from Ref. 11]. 



meaningful, one needs to choose the phase convention for a to be the same as that for a. I 
use a 2 and d 2 to match the phase conventions of the experiment and the theory. I calculate 
d 2 using the same operator used in calculating fio; the d term does not contribute here. 12 
I will fix the phase convention in this way in the next section, followed by the comparison 
between a and a . 



V. RESULT 

A. 1 D 2 —>■ 3 Si transition amplitude 

I calculate the 1 D 2 — > 3 S\ transition amplitude, a^, with the operators presented in 



Sec. mi Then, I compare a.2 with a 2 from the data ll|]. In calculating a 2 , I set the sharp 
cutoff A = oo because there is no counter term at this order which takes care of the high 
momentum components of the operators integrated out. In Fig. |2j I show d 2 obtained with 
the CD-Bonn A^A^ potential, as a function of rj (= q^^/m^). As stated in Sec. 1111} I 
consider neither the Coulomb interaction nor isospin violation effects. I use the CD-Bonn 



12 Because of the finite cutoff, the d term gives a very small contribution, which I safely neglect here. 
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potential for the proton-neutron channel. The solid curve is obtained with Ha = 2.10 (and 
the corresponding q's) while the dashed curve with Ha = 2.68. I choose the overall sign 
for a,2 such that a 2 and d 2 have the same sign, thereby fixing the phase convention. The 
phase convention does not change if I use a finite cutoff such as A = 800 MeV in calculating 
d 2 . As observed in Fig. [2J, d 2 with Ha = 2.10 falls exactly on the experimental data. We 
see the dependence of d 2 on the choice of the NN potential in Fig. [3j I used only Ha = 
2.10 in Fig. [3j Although there is some dependence on the A^A" potential, all a 2 are on the 
vicinity of the data. ;^PT gives a successful description for a 2 at this order. Although a 2 
is significantly larger than a® in magnitude, I stop the discussion on a 2 here and will give 
a detailed discussion on ao in the following subsections. This is because I am primarily 
concerned with the validity of the bridging program in this work, and ao is the one to be 
examined for this purpose. 



B. 1 S'o — ► 3 «Si transition amplitude 

Now I move on to the 1 So — > 3 Si transition amplitude, 5q. Because the relative phase 
between d 2 and ao is fixed within the theory, and the overall phase has been fixed by 
comparing a 2 and a 2 , I am now able to compare ao and do with the same phase convention. 
It is noted that my phase convention for ao is the same as the convention used in Ref. Q|. 

At first, in order to see the importance of the A, I mention a result obtained with the 
7r-production operator without the A. I use the operators given in Eqs. (J5J), ffTUl) . (fT51) - ([TUj) 
with the parameters taken from Ref. [sj . The cutoff function is also the same as that used in 
Ref. [sl (the Gaussian cutoff). I found a very large cutoff dependence of a calculated with 
the A-less x?T. Depending on the cutoff (A x = 500, 600 and 800 MeV), the contributions 
from the two-body operators to a are different by a factor of 4, which demonstrates the 
failure of the A-less in describing the p-wave 7r-pro duct ion. I note that this result is 
quite different from the situation of the low-energy weak process where the A-less %PT gives 
the cross sections with a small cutoff dependence. 

Next I present results obtained with the operator including the A. The operators have 
been presented in Sec. [TT1 and the LEC d has been fixed using the low-energy weak process in 
Sec. IIHI I am interested in how reasonably and reliably one can predict ao for the pp — > pmr + 
reaction using this xPT-based operator. For this purpose, I examine the dependence of ao 
on several inputs, A, 1%a, C 2 A , and the A^A" potential. At first, in Fig. HJ I present a result 
obtained with the CD-Bonn A^A^ potential and two choices of the Ha value. The partial wave 
amplitude So is rather different from the experimental counterpart, ao- For comparison, I 
also show a result obtained without the d term. This result with d = is similar to the case 
with 5 = in Fig. 3 of Ref. \^ where a negative d value (e.g., 5 = —0.2) brings the theoretical 
amplitude into the agreement with the experimental one. In my calculation, however, the 
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FIG. 4: The transition amplitude do for pp — > 
pmr + . The chiral NLO 7r production operator 
and the CD-Bonn NN potential are used. The 
solid and dashed curves correspond to Ha = 2.10 
and 2.68, respectively; A = 800 MeV. The upper 
(lower) curves are obtained with (without) the 
d term. Experimental data are from Ref . . 
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FIG. 5: The transition amplitude do for pp — > 
pmr + . The CD-Bonn NN potential is used. The 
solid, dashed and dotted curves correspond to 
A = 800, 600 and 500 MeV, respectively; h A 
= 2.10. The upper (lower) three curves are ob- 
tained with (without) the d term. Experimental 
data are from Ref. 11]. 



positive d value (see the last row of the second and third columns in Table H]) separates do and 
a>o furthermore; even the sign of d fixed by the low-energy weak process is inconsistent with 
the experimental data of the 7r-production. I change the values of A and Ha and examine the 
dependence of So on these inputs. The cutoff dependence is shown in Fig. The situation 
of the disagreement does not change. 13 This result shows that the bridging program among 
reactions with quite different kinematics is not necessarily successful. This is understandable 
if we recall another case. For example, a chiral nuclear force 26|, |29| describes the elastic NN 



scattering over a fairly wide energy region, partly because the LECs in it have been fixed 
using data from the same energy region. Therefore, it is no wonder to find that the operator 
fixed in the low-energy process cannot describe well the intermediate-energy process. In 
order to accurately describe the two reactions in the different energy regions simultaneously, 
data from both of the energy regions would be necessary to fix the LECs. It is also expected 
that higher order terms are necessary to accurately reproduce the data from the wide energy 
region, as in the case of the nuclear force. 

We find from Fig. El the cutoff dependence (~10% level for ao; the d term included) which 
is much reduced compared with the A-less case. We also find a certain amount of dependence 
on the choice of Ha (Fig.Hj), even if the d value has been adjusted to eliminate the dependence 



13 For some combinations of h A , A and NN potential, the sign of d is consistent with the data. However, 
the strength is not enough to bring ao into agreement with oq. 
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TABLE II: Contribution from each component of the operator to the matrix element for pp — > dir + 
(second row) and pp — ► de~v e (third row). Each contribution is divided by the "sum", and 
therefore, "IB" actually gives the matrix elements with the same sign for both cases. The CD- 
Bonn potential, Ka = 2.10, and A = 800 MeV are used. Some operators, not shown here, contribute 
to the 7r production by a small amount; at most, "WT" contributes by ~ 0.08. 



kinematics 


IB 


A7T 


ACT 


C3 


c 4 


d 


sum (arb. units) 


pp — > dir + 


-0.60 


-0.54 


0.97 


-1.11 


0.71 


1.57 


-0.0251 


pp — ► de~v e 


0.996 


0.009 


-0.008 


0.015 


0.005 


-0.017 


0.757 



on it at the low-energy kinematics. 14 This means that the A operator and the d contact 
operator have quite different dependence on the kinematics, and it is important to take care 
of each component of the operator individually. The dependence on the kinematics is also 
found in Table HT1 where I tabulated contributions from each component of the operator to do 
for two kinematics; one for low-energy weak process (pp — > de~v e ) where d is fixed, and the 
other for the it- production (pp — > dir + , r\ = 0.5). This situation is in contrast to the bridging 
program done in Ref. sj]. In Ref. jsj], the operator fixed by a low-energy weak process (i.e., 
the tritium /?-decay) was used in another weak process which takes place in a relatively 
similar kinematics, and the result was given with a small cutoff dependence. I might say 
that it is not important to take care of each component of the operator individually in this 
case. Rather, one needs to take care of the sum of each component. This is also a reason 
why the pionless theory, whose two-body operator is just a contact operator, can reasonably 
describe several weak processes in a low-energy region [30t] . 

From Table [Til we a ls° nn d that the contributions from the NLO operators is compara- 
ble to those from the LO operator. The LO contributions mostly cancel each other. For 
this reaction, there is no sign of the convergence in the chiral expansion up to this order. 
Furthermore, there is also a significant cancellation among the NLO contributions, which 
may make the amplitude sensitive to the higher order contributions. This situation reminds 
us of the result in Ref. 
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32| where the pp — > ppir reaction near threshold was studied 
with %PT. They found that some higher order two-pion-exchange mechanism contributes 
more than lower order one-pion-exchange mechanism, leading to a poor convergence of the 
chiral expansion. It would be important to do a higher order calculation of the p-wave ir 
production to see the convergence of the chiral expansion. 

We find in Table HU that the contribution from the d term is substantial. This result may 



1 The variation of Ha should be compensated for by the change of Cj and other higher order one-pion 
rescattcring diagrams; not by the change of d. I will come back to this point later. 



19 




3j 

Q 



0.5 



-0.5 



0.1 



0.3 



0.5 



0.7 



FIG. 6: The transition amplitude do for pp — > 
pnn + . The solid, dashed and dotted curves are 
respectively obtained with the CD-Bonn, AV18 
and N 3 LO NN potentials; h A = 2.10 and A = 
800 MeV. The upper (lower) three curves are 
obtained with (without) the d term. Data are 
from Ref. 11]. 



FIG. 7: The transition amplitude do for pp — > 
pmr + . The solid, dashed, dotted and dash- 
dotted curves are respectively obtained with the 
CD-Bonn, AV18, Reid93 and Nijmegen I NN 
potentials; h A = 2.10 and A = 500 MeV. The 
N 3 LO potential almost gives the solid curve. 
The d term is included. Data are from Ref. 111 ]. 



also indicate the importance of a higher order calculation, which I will argue in the following. 
The large contribution from the d term means that its dependence on the kinematics is rather 
influential on ao- Because I am working at the NLO, two-pion-exchange (TPE) mechanisms 
are not explicitly considered but mimicked by the d term. If the TPE and the d term have 
rather different dependence on the kinematics, So given by the NLO and NNLO calculations 
may be significantly different. 

I show results obtained with various NN potentials to see the dependence of a on it. 
The result is shown in Fig. [61 We find a considerable dependence on the NN potential, even 
though all NN potentials give almost the same amplitude in the low-energy region where 
the d value has been fixed. The d term is quite sensitive to the short- distance behavior of 
the wave function because of its point-like nature. Therefore, when the cutoff is relatively 
large and thus the short-distance behavior of the wave function is very dependent on the AW 
potential, the d value is also largely dependent on the AW potential. In case of A = 800 MeV, 
for example, a part of the dependence of ao on the AW potential is ascribable to that the 
d term alone is too simple to compensate for the difference in the short distance physics. 
However, I consider that there is a more influential source of the dependence on the AW 
potential, because we still find the dependence in the case of A = 500 MeV (Fig. EJ). Let me 
explain more on this in the following. In the rest of this paragraph, I discuss the case with 
A = 500 MeV. At this cutoff, the details of the short distance physics have been integrated 
out substantially, and therefore there is no significant difference among the wave functions 
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for different iViV potentials any more. Recall that low-momentum NN potentia 
from various phenomenological iViV potentials are very similar at A ~ 400 MeV 
this situation, the d values for different iViV potentials should be almost the same to give 
almost the same So- Contrary to this expectation, I obtained quite different d as seen in 
the second row of TABLE [J leading to rather different So- 15 I suspect that the d term 
alone is too simple to simulate the operator to be simulated, 16 and that my procedure of 
determining d, discussed in the previous section, gives the d term which is far from being 
as approximate to the operator to be simulated as it can be. I suspect that this is the main 
source of the NN dependence of So- The relatively large /^-dependence of d (and So) is 
also likely to have the same origin. In order to improve the situation, one would need to 
do a higher order calculation, and include a few more contact operators so that a better 
simulation can be done. 17 

Because I treated the contact AW — > A^A interaction in the phenomenological manner, 
as discussed in Sec. HIJ it is informative to study an impact of changing the coupling on the 
amplitude S . For this purpose, I change the C^ A value by ± 10%, re-fit the d value in the 
way discussed in Sec. IIII[ and calculate So- I found that So is less dependent on the variation 
of C^ A than that of h^- This result may be understood as follows. The contact-induced 
A-excitation mechanism [Fig. ED^c)] is similar to the d term [Fig. [T](e)] in the sense that 
Fig. [T](c) is reduced to Fig. [T](e) in the A-less theory. Therefore, one may expect that the 
variation of C^ A is fairly well compensated for by the change of d. 

C. 1 Sq — ► 3 Si transition amplitude with one more counter term 

As seen in the previous section, the bridging program was not successful; the xPT-based 
operator with d(= m N f%d) fixed by the low-energy weak process does not reproduce the 
partial wave amplitude for the 7r-production, ao, extracted from the data. Because I have 
pointed out several reasons for expecting a higher order calculation, here I try to see what 
happens there by doing a simple extension of my calculation. The extension is to add a 

15 In the next subsection where I include one more contact term, we will see that the expectation is realized. 

16 Even though we do not know, there should exist an operator which the d term tries to simulate. The 
operator should include all non-negligible chiral operators other than those explicitly considered already. 
I refer to the operator as the operator to be simulated. 

17 One might suspect that the dependence on the NN potential is due to a difference in the phase shift of 
the partial wave in the energy region above the pion production threshold. Among the NN potentials 
I use, the phase shifts from the Reid93 and the N 3 LO potentials are noticeably different from the others 
at this energy. However, I do not consider the difference in the phase shift to be influential because the 
inclusion of one more contact term results in that all of the NN potentials with A = 500 MeV give 
essentially the same ao over the entire energy region under consideration. Note that the contact term 
does not reflect the difference in the long range behavior (phase shift) of the wave function. 
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TABLE III: Dimensionless contact couplings, d and e. The first column is the sharp cutoff value, 
and the others are values of the LECs. The 2-4th rows are for the d values while 5— 7th rows for 
the e values. For each NN potential, the left side is the LECs for hj\ = 2.10 while the right side 
for h A = 2.68. 

A CD-Bonn AV18 Reid93 Nij I N 3 LO 



500 


0.57 


0.57 


0.45 


0.46 


0.43 


0.44 


0.51 


0.51 


-0.22 


-0.21 


600 


0.49 


0.48 


0.20 


0.19 


0.15 


0.14 


0.37 


0.36 


-4.12 


-4.08 


800 


-1.15 


-1.34 


3.06 


3.30 


3.12 


3.35 


-5.41 


-6.14 


8.42 


8.35 


500 


0.85 


0.85 


0.94 


0.94 


0.96 


0.95 


0.90 


0.90 


1.44 


1.43 


600 


1.10 


1.11 


1.35 


1.36 


1.40 


1.41 


1.20 


1.22 


5.15 


5.12 


800 


3.39 


3.60 


-1.55 


-1.84 


-1.62 


-1.89 


8.38 


9.23 


-7.83 


-7.75 



higher order counter term to the NLO xPT-based operator used in the previous section. I 
use the following counter term of 0[(m n /m^) 2 ] 18 : 

= --JL^rtr -a-VnN (N^ 2 N + h.c.) , (48) 

where the dimensionless LEC is denoted by e. This is not a general form of the counter terms 
at this order, and there are other counter terms with different spin-isospin and derivative 
structures. However, it is sufficient to consider only this counter term for my purpose of 
gaining a rough insight into a higher order calculation. 

Now I have the two independent LECs: d and e. I fix these two LECs so that the 
following two conditions are satisfied. The first condition is the same as that used in Sec. II III 
for fixing d. This condition is from the low-energy weak process. The second condition is 



due to ao extracted from the experimental data for the pp — > pnn + reaction 11] . I choose the 
LECs so that the central value of the empirical amplitude, a^i] = 0.5) = —0.418 fib 1 ^ 2 , is 
reproduced. As seen in the previous section, one cannot reproduce the 7r-production partial 
wave amplitude (a ) by using d fitted to the amplitude of the low-energy weak process alone. 
The disagreement between ao and ao is rather serious; even the sign of d is inconsistent with 
the data in some cases. Therefore, it is not obviously expected that the addition of the e 
term with a natural strength brings ao into agreement with ao- However, such a set of d and 
e does exist, as presented in Table [LTTI 



18 There exist ©[(tott/toat) 3 / 2 ] operators coming from pion-loop diagrams. For a simple analysis, I use the 
0[(m 7r /mN)' 2 } counter term. 
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FIG. 8: The transition amplitude do for pp — > 
pmr + . The chiral NLO 7r production operator 
plus the counter term [Eq. (fl8|) ] and the CD- 
Bonn AW potential are used. The solid, dashed 
and dotted curves correspond to A = 800, 600 
and 500 MeV, respectively, andjiA = 2.10. Ex- 
perimental data are from Ref. 
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FIG. 9: The transition amplitude do for pp — > 
pmr + . The AV18 AW potential is used. The 
other features are the same as Fig. |HJ 



III. 



I calculate a with the parameter sets in Table IHIl The results for the CD-Bonn and 
the AV18 AW potentials are shown in Figs. [8] and [9j respectively. I show results obtained 
with tiA = 2.10 only. The dependence of d on is small in this case. By construction, 
the calculated a goes through the central value of the data at r\ — 0.5. The 77 dependence 
of ao is consistent with the data for all AW potentials and for all A. Particularly, all AW 
potentials give essentially the same do for A = 500 MeV. From Table IIH] we also observe that 
the couplings d and e for A = 500 MeV are similar for all the AW potentials, except for the 
N 3 LO potential. In fact, this result has been expected; see footnote [151 For A = 800 MeV, 
a relatively large dependence of ao (two types of the behavior) on the AW potential is 
observed, within the consistency with the data. Regarding the CD-Bonn and the Nijmegen 
I potentials, a (A = 800 MeV) is very similar to a (A = 500 or 600 MeV); see Fig. E 
On the other hand, the other AW potentials give do (A = 800 MeV) whose 77-dependence 
is quite different from ao (A = 500 or 600 MeV); see Fig. [9j From a viewpoint of the 
renormalization group, the operators with A = 500, 600 and 800 MeV should be related by 
integrating out the high momentum states whose effects are simulated by the renormalization 
of d and e. If this relation is realized, these operators should give the same observables, to 
the extent that the contact operators simulate the high energy modes integrated out. From 
this viewpoint, the results for the CD-Bonn (Fig. [8]) and the Nijmegen I are understandable, 
while the others (Fig. [9]) are not. However, the situation may change in a correct higher 
order calculation where the TPE is explicitly considered. This is because the model space 
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employed here probably has a resolution at which the contact operators cannot accurately 
simulate the intermediate-range mechanism such as TPE. If the TPE is explicitly considered, 
then the contact operators do not have to mimic the intermediate-range mechanism, and 
more accurately describe the shorter-range mechanism. 

I look into the contact operators for different NN potentials, and understand the similar 
(different) ^-dependence of do for A = 500 (800) MeV among different NN potentials. As 
representatives, I show the contact operators for the CD-Bonn and the AV18 potentials in 
Figs. [TD] and dU respectively. We can find that the contact operators, including both d and 
e, for A = 500 (800) MeV are very similar (different), leading to the quite similar (different) 
rj dependence of clq. Starting with the operators for A = 800 MeV, one can reduce A to 
examine the running of the operators using the Wilsonian RG equation. I refer the readers 
to Refs. [35j, |36| for a detailed explanation of how I reduce A, and just show the result here. 



In Fig. [T21 the RG running of the contact operator for the CD-Bonn potential is given. For a 
comparison, the contact operator for A = 500 MeV, whose LECs are given in Table 1111} are 
also shown. We can see that the contact operators (A = 500 MeV) derived in the different 
two ways, one from directly fitting to the data and the other from the RG running, are fairly 
similar. Although the figure shows the result for the diagonal components of the operators, 
the same level of the agreement is confirmed for off-diagonal components. For the AV18 
potential, however, the two contact operators (A = 500 MeV) with the different origins are 
quite different as shown in Fig. [131 This result has been expected by observing Fig. [9] where 
the contact operators with A = 500 and 800 MeV do not give the same result, indicating 
that the two operators are not equivalent. 

I compare the contact operator including both d and e (Table UTTj) . and those with only 
d (Table [I]) in Figs. [10] and [TTJ A naive expectation is that a contact operator with d only 
(a dashed curve) is an approximation of the one with d and e (the solid curve). In the two 
figures, this expectation is not always the case, and therefore we have to be careful about 
the convergence of the chiral expansion of the short distance physics. In any case, the d 
term alone cannot be a good approximation of the d plus e terms and, for that matter, not 
a good approximation of the operator to be simulated. It is recalled that this observation 
has been used in the previous subsection to argue over a reason for the dependence of a 
on Ha and the NN potential. By including the e term, the situation is much improved in 
this point. We found the very small dependence of do on the choice of Ha and on the NN 
potential for A = 500 MeV. Regarding the d and e values, as seen in Table IIII1 they hardly 
depend on Ha- This is quite consistent with the discussion given in footnoted] 
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FIG. 10: The contact operators for the CD- 
Bonn potential. The diagonal components in 
the momentum space are given. The solid curves 
are parameterized by both d and e, while the 
dashed curves by d only. The ^-coordinate value 
at the end point of a curve indicates the value 
of A for the operator. Ka = 2.10. 



FIG. 11: The contact operators for the AV18 
potential. The other features are the same as 
Fig. DIH The scale is also the same as Fig. [TOl 



VI. SUMMARY 



I determined the LEC d using the low-energy weak process and then used it to predict 
the partial wave amplitude, ao ( 1 5o — > 3 5*i), for the pp — > pnn + reaction. Through this 
work, I tried to explore the power of %PT that enables one to bridge different reactions. 
My investigation is more stringent test of this aspect of %PT than similar analyses in the 
literature because the reactions connected through xPT here take place under significantly 
different kinematics. 

I started with the chiral Lagrangian including the nucleon, pion and A. It is mandatory 
to include the A explicitly for describing the p-wave ir production. With the Lagrangian, I 
constructed, up to NLO of the chiral expansion proposed in Ref. (5), [ljj, a set of operators 
which describes the Gamow- Teller transition in low-energy weak processes and the p-wave 
7r-production. I fixed the remaining unknown LEC d (indirectly) using the experimental 
data of the low-energy weak process. Then I calculated the partial wave amplitudes (a) for 
the pp — > pnn + reaction. I chose the phase convention such that a<i (}D2 — ► 3 S\) has the 



same phase as 02 extracted from the data in Ref. llj. My prediction of ao using the NLO 
operator does not agree with the data. I used several different sets of the inputs such as the 
AW potential, Ha and A. Even though there is some dependence of ao on these inputs, all 
results differ from the data in the similar manner. Unfortunately, even the sign of d, fixed by 
the weak process, is sometimes not consistent with the data. This result indicates that the 
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FIG. 12: RG running of contact operators for 
the CD-Bonn potential. The solid curves are 
the same as those in Fig. [H)J (The curve for 
A = 600 MeV is not shown.) The dashed and 
dotted curves are derived from the solid curve 
for A = 800 MeV by solving the RG equation. 
The dashed (dotted) curve is obtained with the 
kinematics for the pp — > dir + (pp — ► de + u e ) re- 
action. 
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FIG. 13: RG running of contact operators for 
the AV18 potential. The solid curves are the 
same as those in Fig. [TTJ The other features are 
the same as Fig. [T2J 



bridging program between the two reactions with quite different kinematics is not always 
successful. In the literature, we have sometimes found an argument which supposes that the 
bridging program works. Given the result here, it is clear that we need to seriously study 
a feasibility of the bridging program, particularly for reactions with different kinematics. 
This conclusion may be disappointing, but still understandable if we recall the success of 
the chiral nuclear force. The chiral nuclear force accurately describes the NN scattering 
over a wide energy region, partly because the LECs included have been fixed using data 
from the same energy region. In fact, there are several reasons to expect a higher order 
calculation to significantly improve the situation. First of all, one may naively think that 
the NLO, tree level, operator is too simple to describe the 7r production. We know that 
two-pion-exchange mechanism and higher order contact terms are necessary for accurately 
describing the AW elastic scattering near the ir production threshold. Second, the d term 
largely contributes to So, and therefore it would be important to extract the TPE mechanism 
from the d term, thereby describing the intermediate-range mechanism more elaborately and 
reducing the role of the d term. Third, a is considerably dependent on the choice of the 
NN potential, which means that the single contact term is too simple to cancel out different 
short distance physics for different NN potentials. I also argued that a substantial part of 
the W-dependence (and /i^-dependence) of do is likely to be ascribable to the fact that 
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the d term alone is too simple to simulate the operator to be simulated. Meanwhile, a higher 
order calculation is also desirable to see the convergence of the chiral expansion; regarding 
the pp — ► pnir + reaction, we found no sign of the convergence up to this order. 

In order to explore, even roughly, a result of a higher order calculation, I added a higher 
order counter term, with the LEC e, to the NLO operator. I fitted the LECs d and e to both 
the Gamow- Teller amplitude for the low-energy weak process and do for the tt production. I 
found a set of the LECs with the natural strength. The LECs are mostly independent of the 
choice of Ha as they should; without the e term, however, d is rather dependent on Ha- With 
this parameter set, //-dependence of do is described in a way consistent with experimental 
data, irrespective of the choices of the NN potential, Ua and A. This result would be an 
indication that a higher order calculation is promising. I found that the single d term is not 
always a good approximation of the d plus e terms. I also showed using the RG analysis 
that the contact operators with different cutoff are not always equivalent. These findings 
tell us to be careful about the convergence of the chiral expansion, and also the importance 
of going to a higher order calculation. 
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APPENDIX A: MULTIPOLE EXPANSION OF A n OPERATOR 

I present our calculational procedure for the AlT operator [Eq. (J5J]. I start with the 
Fourier transform of O^: 

dk 3 _ ifc . r k <t 2 ■ k /(A|) 



(2tt) 3 m! 2 + k 2 m A - m N - p 2 /m N + (p' + q 7r /2) 2 /2/i 

where I only consider one term in Oa-k', the other terms and a constant factor are omitted. 
The function f(k) is a cutoff function: f(k) = exp (— k 2 / K 2 G ) with Aq = 2 GeV. In this 
equation, I expand the energy denominator as follows: 



m A -m N - p 2 /m N + (p' + q 7T /2) 2 /2^ 
4nY / (-lYV2lTT[Y i (p') ® Yeiq^Dtij/,^) , (A2) 
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where Dp{p' , q n ) is defined by 

Dt(p',q*) = - J_ D{p ,q^)P E (z)dz 



2 J-i 
1 f l P t {z)dz 



2 J -i (j/qj2fi)(J3 + z 
(-1) 



ip'q^/2/i) 



Qtifl) , (A3) 



with 



= m A ~ m N -p 2 Jm N + (p' 2 + g|/4)/2/i _ 
p'q n /2fi 

and z = p' ■ q n . The function Qi{j3) is the Legendre function of the second kind, and is given 
by 

I write Eq. (jSH) using the expanded form Eq. flA2j) . When I retain only the first term (£ = 0) 
of the expansion, I obtain 



(27r) 3 + /c 2 p'qv 

i^o(r) + ^ C^( # ) ® ^(D 4^^o(/3) , (A6) 

O O / P Qtt 



with 



1 r 00 A* 
F (r) = — / dk—— —Ukr)f(k) , 
2ir z Jo k + m[£ 

1 r°° k 4 

I take a matrix element of Eq. ( 1A1I) after setting q n = 0, and compared it with the matrix 
element of Eq. (1A6|) in which q n (^ 0) is fixed by the kinematics. In the kinematical region 
of interest, I found a small correction (~ 1.5%). The use of Eq. ( 1A6I) may be regarded as 
an inclusion of the lowest order correction [0(g 2 )] from finite q n , even though there is still 
another 0{q 2 ) correction. I do not consider the higher order {£ > 1) terms in Eq. (1A2I) to be 
influential on our result because: the 0(q 2 ) correction from Eq. (1A6I) is small; the expansion 
in Eq. (]A2j) may be regarded as an expansion in terms of z//3, and z//3 -C 1 in most of the 
kinematical region of interest. I use Eq. (1A6I) in our calculation. 
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